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DISCOVERY OF HARD X-RAY DELAYS IN THE X-RAY EMISSION OF THE SEYFERT 1 GALAXY MRK 110: 

POSSIBLE EVIDENCE FOR COMPTONIZATION 
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ABSTRACT 

We report the discovery of hard X-ray delays in the X-ray emission of the Seyfert 1 galaxy Mrk 110, 
based on a long XMM-Newton observation. Cross correlation between the X-ray light curves of 
different energy bands reveals an energy dependent delay ranging from a few minutes to an hour. 
We find that the energy spectrum can be modeled by Comptonization of disk blackbody photons. 
The energy dependent delay can be modeled as due to the effect of Comptonization in a hot plasma 
confined within 10 Schwarzschild radius of the black hole. We discuss our results in the context of 
inverse Comptonization of the soft photons by highly energetic plasma. 

Subject headings: galaxies: active — galaxies: individual (Mrk 110) — radiation mechanisms: general 
— accretion, accretion discs — X-rays: galaxies 



1. INTRODUCTION 

The dominant X-ray radiation mechanism in accreting 
black holes is commonly thought to be inverse Comp- 
ton scattering of low energy photons by a cloud of hot 

I T ^ 10^ — 10^ K) electron plasma near the black hole 
Sunvaev & TitarchuklllQSdl) . The physical geometry of 
the coronal plasma responsible for scattering the pho- 
tons is not well constrained by the current observations. 
The ubiquity of a disk blackbody component accompa- 
nied by a power-law tail in the overall spectra of galactic 
black hole candidates (GBHCs) motivated several work- 
ers to develop the so c alled two phase accretion disk- 
corona models (see eef..lHaardt f: Maraschilll9911 1199.'^ 
iPoutanen. Krolik^ fc R ydd Il997t iSvensson fc Zdziarskil 
E994; Ste rn et all 119951 iBeloborodovl \199^ . In these 
models, the blackbody radiation from the cold disk en- 
ter the hot corona and is Comptonized into X-rays. A 
part of the hard X-rays from the corona, being repro- 
cessed in the disk, produce the reflection hump. The 
geometry of the corona controls this feedback mech- 
anism which in turn determines the spectral slope of 
the escaping radiation. The Kompaneets y-parameter 
(and hence the temperature) is determined by the en- 
ergy balance between heating and cooling mechanism 
inside the plasma. One important question in all such 
models is the method by which the gravitational en- 
ergy is converted to the energy of electrons. The 
ideas explored in the literature include magnetic flares 
iHaardt. Maraschi fc Ghisellinil 119941: iStern etld] 119951 
IPoutanen & SvenssonI 119961: IBeloborodovl Il999albl ). ad- 
vection do minated disk very close to the black hole 
ijNaravan fc Yi 1 994.), Bondi type free fall beyond a 
shock region ijChakrabarti fc TitarchukI Il995j) . but, as 
yet, there is no consensus on the exact mechanism. A 
general prediction of the above models is that the hard 
X-ray variations should lag behind those in softer bands, 
as hard photons undergo more number of scatterings in 
the plasma before escape. Measurement of such time-lags 
has the potential of constraining the size of the region 

^ Tata Institute of Fundamental Research, Mumbai- 
400005, India, email: suraiit@tifr.res.in, arrao@tifr.res.in, 
url; |http: / / www.tifr.res.in/^surajit, 



which in turn will help us to understand the mechanism 
by which the electron cloud is energized. 

Frequency dependent time lags have been observed in 
Gala ctic black hole candidates l)Mivamoto et al.l 119881 
1991 as well as in a few activ e gala ctic nuc lei (A GNs; 
Papadakis. Nandra. fc Kazana£ll200ll lZhang|l2002|) . In- 
terpreting these lags as due to Comptonization requires 
the size of the emitting region to be very large, typically 
several thousand Schwarzschild radii (compatible to the 
lowest frequency at which the lag is observed), leading 
to the problem of heating the electron cloud at large dis- 
tances from the black hole. Hence, sometimes these lags 
are interpreted as due t o the energy dep endent asym- 
metries in random shots ijMivamoto fc Kita mato 1989) ; 
that is interpreting the lags as due to the production of 
the variability itself. 

One of the problems of detecting lags at higher frequen- 
cies and relating them to the Comptonization at the in- 
nermost regions of the accretion disk around black holes 
could be observational. If the Comptonization process 
occurs at 10 - 20 Schwarzschild radii and if the Comp- 
tonization process gives a lag which is a factor of few 
larger than the light travel time in these regions, one ex- 
pects a lag of a few tens of milliseconds in a black hole 
of mass 10 M© and about a day in a AGN of mass 10^ 
Mq. Detection of such delays is observationally a diffi- 
cult task. Bright nearby AGNs with a black hole mass of 
10'' M0, where one expects a delay of about an hour, are 
the ideal sources to look for delays due to the process of 
Comptonization. With this motivation, we have searched 
for delays in one of the bright low mass AGN Mrk 110, 
based on a long observation using the XMM-Newton ob- 
servatory. In this Letter we present the cross-correlation 
analysis in different energy bands of this source and show 
that hard X-rays are delayed by a significant amount of 
time (from few minutes to an hour). The 0.3 - 12 keV 
X-ray spectrum can be represented by a Comptonization 
model which can explain the measured hard X-ray delay. 
2. OBSERVATION 

Mrk 110 is a nearby optically bright, radio interme- 
diate {R ~ 1.6) narrow- line Seyfert 1 galaxy (NLSls) 
at a redshift z=0.036. The optical continuum and the 
broad emission lines of Mrk 110 are highly variable 
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Fig. 1.— EPIC-PN light curves (1000 s bin) of Mrk 110 in dif- 
ferent energy range. Start of variability is indicated by downward 
arrows (see text) 



(by a factor of 2 to 8 within a timescale of 10 years - 
IKollatschnv fc Welshll20?nt IKollatschnvll2nO,H) . Mrk 110 

was observed on 2004 November 15 by XMM-Newton for 
47.4 ks. The EPIC-PN cameras were operated in Prime- 
Small- Window observing mode using the Thinl filters. 
Here we use data only from EPIC-PN cameras due to 
the better efficiency, better calibration at lower energies 
and absence of pile-up. Source spectra and light curves 
were extracted from the EPIC images using a circular 
source region centered on the observed source position. 
Background spectra and light curves were derived from 
adjacent 'blank sky' regions. The EPIC spectra were 
binned to give a minimum of 50 counts per bin. The 
XSPEC vll.O and xronos packages were used for spec- 
tral and timing analysis respectively. Errors on fitted 
parameters are quoted at the nominal 90% confidence 
level (Ax^ = 2.7) unless otherwise stated. 

3. TIMING ANALYSIS 

In Figure m we plot the binned light curves of Mrk 110 
in different energy ranges (not corrected for 71% duty cy- 
cle of PN SW mode). The 0.2 - 12 keV light curve shows 
peak to peak variation of approximately 10% within 3 
hours. To quantify the source variability in different en- 
ergy bands, we calculated the fractional variability in 
seven energy bands: Ei (0.2 - 0.3 keV), E2 (0.3 - 0.42 
keV), E3 (0.42 ~ 0.58 keV), E4 (0.58 - 0.8 keV), E5 (0.8 
- 1.2 keV), Eg (1.2 - 2 keV), and E7 (2 - 12 keV), re- 
spectively. The energy bands are chosen such that the 
mean count rate in those bands are approximately same 
(within 10% of average). We find significant variability 
(rms value 2.5-4.5%) in all the energy bands for 500 s bin- 
ning (with a typical error of 0.4%). There is a marginal 
evidence for increasing variability with increasing energy: 
the variability is (2.7 ± 0.4)% for energy <0.6 keV and it 
is (3.6 ± 0.4) % for energy > 0.6 keV. A Structure Func- 
tion analysis l)R,utmanl 119781) of the light curves shows 
that the shortest correlation timescale is more than a 
few thousand seconds. 

A visual examination of the light curves reveals a grad- 
ual decrease in count rate (^ 0.2 cts in 2 hrs) in Ei at 
time 10 ks from the starting time. Similar kind of vari- 
ation is seen after more than an hour in E7. (marked 




-10000 -5000 



5000 10000 
Time Delay (s) 



15000 



Fig. 2.— The CCFs between light curves of 0.2 - 0.3 koV and 
higher energy bands (increasing from top to bottom). The CCF 
values are vertically shifted for clarity, (see text) 



by arrows in figure^. To search for time lags, we cal- 
culated the cross correlation function (CCF) between Ei 
and other energy bands using the crosscor package in 
FTOOLS. The CCFs between Ei and E^ (where i=2,3,..,7) 
are plotted as a function of delay, the successive plots are 
vertically shifted by 1.6, 1.2, 0.9, 0.6, 0.3, respectively 
for clarity. The errors in the CCF values are the stan- 
dard one sigma values which includes only the counting 
statistics errors. The harder light curves systematically 
lag the softer band. To estimate amount of the delay, 
we fitted the central part of the CCF distribution with a 
Gaussian function and derive values of delays for the dif- 
ferent energy bands of 200±700 s, 400±900 s, 900±800 
s, 1800±900 s, 2400±900 s, 4500±900 s, respectively. 
The errors in the delays are estimated by the fit- 
ting method with the prescription A^^ — 4 (for three 
parameters) . In Fig 13 we plot the derived delay as a 
function of the energy of the hard band, along with a 
Comptonization model (see section 5) fitted to the data 
points (x^/dof ~ 1/4). For the hypothesis that there 
is no delay we get a value for x^/dof of 38/6 and for 
the hypothesis of energy independent delay gives a value 
for x^/dof of 18/5. Hence we can conclude that a de- 
lay is detected and it is energy dependent at more than 
99% confidence level. We must, however, caution that 
the above results are based on statistical considerations 
only and do not include systematics like the shape of the 
variation of CCF with delay. We have also derived cross 
correlation using other combinations (CCFs of Ei, w.r.t. 
E4, where 1=1,2,3,5,6,7) to check whether the results are 
artifact of uncertainty of instrumental calibration of en- 
ergy below 0.3 keV and we find similar results. 

4. SPECTRAL ANALYSIS 

We first modelled the X-ray spectrum of Mrk 110 
with a absorbed power-law with Galactic column den- 
sity Nh = 1.6 X 10^". This provides a fairly good fit to 
the PN data in the 2-12 keV range {x'^/v = 1037/1022) 
giving a photon index F = 1.75 ± 0.01. Inclusion of a 
narrow Gaussian emission line improves the fit signif- 
icantly {x? /v — 997/1020), giving the energy of the 
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E (keV) FiG- 4.— XMM-Newton (EPIC-PN) spectrum in the energy 

range 0.3 - 12 keV and the best fit hybrid thermal/non-thermal 
Fig. 3.— Energy dependence of lag between the Ught curves of Comptonization model, 

energy 0.2 - 0.3 keV and higher energies. 



line 6.42 ± 0.02 keV. Extrapolating the best fit model 
to 0.3 keV shows a huge soft excess in the spectrum 
[x^/i, = 94471/1363). Refitting still provides a rather 
poor fit {'X^ /v — 7350/1363). A power-law plus a black- 
body (generally the soft excess component in AGN is 
modeled by thermal black-body but the origin is still un- 
known) to model the soft excess improves the quality of 
fit but ft is stiU not acceptable {x^/v = 2109/1361). The 
temperature of the blackbody becomes kT = 100 ± 2 
eV. Fitting the data with more realistic thermal ac cre- 
tion di sc spectrum like diskb b ijMitsuda et alJlT98^ or 
diskpn (>Gierli nski et al.lll999t) provides a poor fit to data 
(x^/iy = 1780/1361 and xVj^ = 1782/1360 respectively) 
and the temperature becomes very high [kT — 155±2 eV 
and 159±3 eV respectively). These results would seem to 
reject an origin for the soft excess in terms of unmodified 
thermal blackbody emission. Addition of a power-law in- 
stead of blackbody improves the fit {x^ jv = 1680/1361). 
The values of the power-law indices are 2.47 ± 0.02 and 
1.21±0.04. A model in 0.3 - 12 keV range which consists 
of a broken power-law and 6.4 keV Gaussian emission 
Hue improves the fit slightly {x^jv = 1597/1361). The 
power-law indices are 2.29 ± 0.01 and 1.78 ± 0.01, and 
the break energy is 1.66 ± 0.04 keV. 

Reflection off the surface could also produce strong soft 
excess at low energies. To test this the ioni zed reflection 
model pexriv l)Magdziarz & ZdziarskillT99l was fltted to 
the 0.3 - 12 keV spectrum resulting in a poor fit to the 
data {x^ /dof ^ 2185/1361). The best fitting parameters 
were T = 2.26 ± 0.01, R = 8.14 ± 0.12 and ^ < 10"^. 

The energy dependent time lags arc strongly sugges- 
tive of Comptonization of low energy seed photons by 
a population of high temperature electrons (section 1). 
The COMPTT code ((Titarchuklll994tl was used to model 
Comptonization of soft photons in a thermal plasma. A 
power-law plus a Comptonization component give a good 
fit to the data {x^/v = 1477/1359) with T = 1.51 ± 0.05 
and the seed photon temperature kT^b = 69 ± 2 eV. 
The temperature and optical depth of the Comptonizing 
plasma are strongly covariant parameters, and thus can- 
not be constrained si multaneously. Fitting the da ta with 
the compps model llPoutanen fc Svensso nl lT99^ along 
with the power-law gives a reasonable flt {x^/dof = 
1483/1359). In this model the soft flux is equal to the 
sum of the absorbed incident flux from the corona and 
the flux due to local energy dissipation in the cold disk. 



The spectral shape of the soft components are assumed 
to be Planckian, with temperature Tbi, and Tdisk , respec- 
tively {Thh > Tdisk)- The inner disc temperature kTdisk is 
fixed at 40 eV (calculated for black hole of mass 10^ Mq 
assuming standard accretion disk). The best fit values 
of the parameters are r = 3.2lJ;y, kTe = 14.0tf;2 keV, 

fcTfefa = 991^ eV, F = 1.57t°;°t 

The Comptonization model along with a power-law de- 
scribed above has the two spectrally identified continua 
originating in two distinct thermal Comptonizing plas- 
mas. An alternative is that the whole spectrum is pro- 
duced by a single plasma wi th a hy brid thermal/non- 
thermal electron distribution ljCoppil ll999). The definite 
detection of a non-thermal Comptonization component 
re quires high energy and hig h quality data. As noted 
by 'Pounds. Do ne fc Osborng (|1995|) , the X-ray spectra 
of ultrasoft Seyfert galaxies do resemble that of Cyg X- 
1 in its high/soft state. Thus it seems reasonable to 
test whether a hybrid thermal/non-thermal plasma is 
a viable model for the X-ray continuum for Mrk 110. 
To test this idea we tried the hybrid Comptonization 
mo del compps (for a det ailed description of parameters 
see IZdziarskiet aLll200l . We assume a spherical ge- 
ometry of hot plasma. We fltted the model [x^ /dof = 
1455/1359) and the best fit values of the parameters are 
r = 4.8^2:1, kTe = ll-l-?;8 keV, kT^b = 105^3 eV, 
F„, = 2.43t°;;^, 7^„, = 1.19l°:J|. The unfolded spec- 
trum is shown in Figure 0] Similar attempt has been 
done to model the X-ray s pectrum of NLSl galaxy Ton 
S180 l|Va,ugha,n et a,ljl20f)l . 

5. DISCUSSION AND CONCLUSION 

We find significant energy dependent delay between the 
hard and soft X-ray emission in the Seyfert 1 galaxy Mrk 
110. Similar delays were also found in the br ight Seyfert 
1 gala xy MCG-6-30-15 llPonti et al.ll200l . IGallo et alJ 
(j2Q^) found alternating leads and lags in the NLSl 
IRAS 13224-3809. Though other interpretations like ge- 
ometric effects and energy dependent shape of shots can 
also explain the observed energy dependent delays in Mrk 
110, we try to interpret the results as due to Comp- 
tonization, particularly because of the fact that the spec- 
tral analysis favours a two component Comptonization 
model. We consider a static Compton cloud with opti- 
cal depth tt, and electron temperature Q — kTf.jm^.c? . 
A soft seed photon of energy Eq injected into the cloud 
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increases its energy by a factor oi A — 1 + 40 + 160^, 
on an average after each scattering, so that after n scat- 
tering its energy is En = A"Eq. The photon mean free 
path is A « R/max{l,TT) (where R is the size of the 
X-ray emitting region). The time difference between 
successive scatterings is then tc — (i?/c)/max(l, tt), 
where tt = NarR is the Thompson optical depth, and 
(Tt is the cross section of Thompson scattering, N is 
the electron number density of the scattering medium. 
The time n eeded to reach the energy i?„ is tn = ntc 
(jSunvaev fc Tita rchuk 1980; Payne 1980J. We have cal- 
culated time lags of different high energy bands with re- 
spect to 0.2 - 0.3 keV band. Hence in our case Eq = 0.25 
keV. We fit the above equation {tn = ntc) to the result 
(hard lag as a function of energy) of our analysis and 
find that the data is in good agreement with the equation 
(Fig. 01 • Using the values of the best fit parameters of the 
above result and the parameters of the Comptonization 
model fitted to spectrum we get R ~ 2x10^^ m fO i?5. 
This value of R is physically realistic as most of the en- 
ergy is dissipated within f Rs- Bott cher fc Lia ng (f 998} 
pointed out that any scenario in which the observed hard 
time lags are purely due to static Comptonization re- 
quires that the radial extent of the hot corona exceeds 
~ 10^ Rs of a solar mass black hole {R > 10^ m). This 
is incompatible with current models of accr etion flow s 
onto galactic black holes (see the review byEang 'l998) 
even from simple energy arguments (see section 1). But 
the results reported here pertains to lags at very short 
times scales and hence is consistent with the static Comp- 
tonization model. The harder power-law emission ex- 
tending to 12 keV (and presumably beyond) can also be 



produced by Comptonization, either in another purely 
thermal plasma or non-thermal electrons in a plasma 
with a hybrid thermal/non-thermal distribution. 

The origin of the hot plasma above the accretion disk 
is quite debatable. If the radiation pressure inside the 
accretion disk very close to the black hole is very high 
then the accretion disk can be puffed up and the hot 
thermal plasma can be formed just above the accre- 
tion disk. Hybrid thermal/non-thermal plasmas have of- 
ten been successfully u sed to model the observed data 
((Cicrhnski et al. 1999: Poutanen fc Coppi.199 8). A pos- 
sible origin of the non-thermal component is in pro- 
cess of ma gnetohvdrodvnarn ic turbulence occurring in 
the corona lILi fc Milleilfl997l) . Stochastic acceleration of 
particles by stochastic gyroresonent acceleration in ac- 
creting plasma can also accelerate p articles to higher en- 
ergies ( Dcrme r^etallll996t lETet al.l ll996') . 

We conclude that inverse Compton scattering of soft 
photons by highly energetic electron distribution pro- 
vides a satisfactory explanation of the hard X-ray time 
lag observed in the XMM-Newton observation. The en- 
ergy spectrum in the 0.3 - 12 keV band can be modelled 
either by two component Comptonization or by a hybrid 
thermal/non-thermal Comptonization. It is not possible 
to distinguish between them without good quality high 
energy data. 

This work is based on observations obtained with the 
XMM-Newton, an ESA science mission with instruments 
and contributions directly funded by ESA member states 
and the USA (NASA). Authors are grateful to H. Netzer, 
G. Dewangan and S. Mandal for useful suggestions. We 
thank the referee for useful comments. 
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